Survey of mitochondrial sequences integrated into the bovine nuclear genome.

Nuclear copies of the mitochondrial DNA (NUMTs) have already been described in several species. In this context, we identified and analysed 166 bovine NUMT regions with a total length of 430 kbp, representing about 0.02% of the cattle nuclear genome. Copies of all mitochondrial regions were detected in the nuclear genome, with distinct degrees of sequence similarity to the mitogenome. Some NUMT regions include large mitogenome segments and show high similarity to the organelle genome sequence. NUMT regions are frequently modified by insertion of repetitive sequences and by sequence rearrangements. We confirmed the existence of 29 NUMT regions by PCR amplification using DNA from the cow (Dominette) which was used to generate the bovine genome reference sequence, ruling out the possibility that these NUMTs could be artifacts of the genome assembly. As there are NUMT regions with high similarity to the mitogenome, special care is needed when designing primers for mitochondrial DNA amplification. Our results can therefore be used to avoid co-amplification of bovine nuclear sequences similar to mitochondrial DNA.

Α de novo 3.8-Mb inversion affecting the EDA and XIST genes in a heterozygous female calf with generalized hypohidrotic ectodermal dysplasia.

BACKGROUND: In mammals, hypohidrotic ectodermal dysplasia (HED) is a genetic disorder that is characterized by sparse hair, tooth abnormalities, and defects in cutaneous glands. Only four genes, EDA, EDAR, EDARADD and WNT10A account for more than 90% of HED cases, and EDA, on chromosome X, is involved in 50% of the cases. In this study, we explored an isolated case of a female Holstein calf with symptoms similar to HED. RESULTS: Clinical examination confirmed the diagnosis. The affected female showed homogeneous hypotrichosis and oligodontia as previously observed in bovine EDAR homozygous and EDA hemizygous mutants. Under light microscopy, the hair follicles were thinner and located higher in the dermis of the frontal skin in the affected animal than in the control. Moreover, the affected animal showed a five-fold increase in the number of hair follicles and a four-fold decrease in the diameter of the pilary canals. Pedigree analysis revealed that the coefficient of inbreeding of the affected calf (4.58%) was not higher than the average population inbreeding coefficient (4.59%). This animal had ten ancestors in its paternal and maternal lineages. By estimating the number of affected cases that would be expected if any of these common ancestors carried a recessive mutation, we concluded that, if they existed, other cases of HED should have been reported in France, which is not the case. Therefore, we assumed that the causal mutation was dominant and de novo. By analyzing whole-genome sequencing data, we identified a large chromosomal inversion with breakpoints located in the first introns of the EDA and XIST genes. Genotyping by PCR-electrophoresis the case and its parents allowed us to demonstrate the de novo origin of this inversion. Finally, using various sources of information we present a body of evidence that supports the hypothesis that this mutation is responsible for a skewed inactivation of X, and that only the normal X can be inactivated. CONCLUSIONS: In this article, we report a unique case of X-linked HED affected Holstein female calf with an assumed full inactivation of the normal X-chromosome, thus leading to a severe phenotype similar to that of hemizygous males.

A First Step Toward Unraveling the Energy Metabolism in Endurance Horses: Comparison of Plasma Nuclear Magnetic Resonance Metabolomic Profiles Before and After Different Endurance Race Distances.

Endurance racing places high demands on energy metabolism pathways. Metabolomics can be used to investigate biochemical responses to endurance exercise in humans, laboratory animals, and horses. Although endurance horses have previously been assessed in the field (i.e., during races) using broad-window Nuclear Magnetic Resonance metabolomics, these studies included several different race locations, race distances, age classes, and race statuses (finisher or elimination). The present NMR metabolomics study focused on 40 endurance horses racing in three race categories over 90, 120, or 160 km. The three races took place in the same location. Given that energy metabolism is closely related to exercise intensity and duration (and therefore distance covered), the study's objective was to determine whether the metabolic pathways recruited during the race varied as a function of the total ride distance. For each horse, a plasma sample was collected the day before the race, and another was collected at the end of the race. Sixteen, 15, and 9 horses raced over 90, 120, and 160 km, respectively. Proton NMR spectra (500 MHz) were acquired for these 80 plasma samples. After processing, the spectra were divided into bins representing the NMR variables and then classified using orthogonal projection on latent structure models supervised by the sampling time (pre- or post-race) or the distance covered. The models revealed that the post-race metabolomic profiles are associated to the total ride distance groups. By combining biochemical assay results and NMR data in multiblock models, we further showed that enzymatic activities and metabolites are significantly associated to the race category. In the highest race category (160 km), there appears to be a metabolic switch from carbohydrate consumption to lipid consumption in order to maintain glycaemia. Furthermore, signs of protein breakdown were more apparent in the longest race category. The metabolic shift seen in the different racing categories could be related to a mixture of three important factors that are the ride distance, the training status and the inherited endurance capacity of the various horses competing.

Survey of allele specific expression in bovine muscle.

Allelic imbalance is a common phenomenon in mammals that plays an important role in gene regulation. An Allele Specific Expression (ASE) approach can be used to detect variants with a cis-regulatory effect on gene expression. In cattle, this type of study has only been done once in Holstein. In our study we performed a genome-wide analysis of ASE in 19 Limousine muscle samples. We identified 5,658 ASE SNPs (Single Nucleotide Polymorphisms showing allele specific expression) in 13% of genes with detectable expression in the Longissimus thoraci muscle. Interestingly we found allelic imbalance in AOX1, PALLD and CAST genes. We also found 2,107 ASE SNPs located within genomic regions associated with meat or carcass traits. In order to identify causative cis-regulatory variants explaining ASE we searched for SNPs altering binding sites of transcription factors or microRNAs. We identified one SNP in the 3'UTR region of PRNP that could be a causal regulatory variant modifying binding sites of several miRNAs. We showed that ASE is frequent within our muscle samples. Our data could be used to elucidate the molecular mechanisms underlying gene expression imbalance.

Linkage disequilibrium and past effective population size in native Tunisian cattle.

To carry out effective genome-wide association studies, information about linkage disequilibrium (LD) is essential. Here, we used medium-density SNP chips to provide estimates of LD in native Tunisian cattle. The two measures of LD that were used, mean r2 and D', decreased from 0.26 to 0.05 and from 0.73 to 0.40, respectively, when the distance between markers increased from less than 20 Kb to 200 Kb. The decay in LD over physical distance occurred at a faster rate than that reported for European and other indigenous breeds, and reached background levels at less than 500 Kb distance. This is consistent with the absence of strong selective pressure within the Tunisian population and suggests that, in order to be effective, any potential genome-wide association mapping studies will need to use chips with higher marker density. An analysis of effective population size (Ne) based on LD data showed a decline in past Ne, with a sudden drop starting about eight generations ago. This finding, combined with the high levels of recent inbreeding revealed by runs of homozygosity (ROH) analysis, indicate that this population is endangered and may be in urgent need of a conservation plan that includes a well-designed genetic management program.

An initiator codon mutation in SDE2 causes recessive embryonic lethality in Holstein cattle.

Researching depletions in homozygous genotypes for specific haplotypes among the large cohorts of animals genotyped for genomic selection is a very efficient strategy to map recessive lethal mutations. In this study, by analyzing real or imputed Illumina BovineSNP50 (Illumina Inc., San Diego, CA) genotypes from more than 250,000 Holstein animals, we identified a new locus called HH6 showing significant negative effects on conception rate and nonreturn rate at 56 d in at-risk versus control mating. We fine-mapped this locus in a 1.1-Mb interval and analyzed genome sequence data from 12 carrier and 284 noncarrier Holstein bulls. We report the identification of a strong candidate mutation in the gene encoding SDE2 telomere maintenance homolog (SDE2), a protein essential for genomic stability in eukaryotes. This A-to-G transition changes the initiator ATG (methionine) codon to ACG because the gene is transcribed on the reverse strand. Using RNA sequencing and quantitative reverse-transcription PCR, we demonstrated that this mutation does not significantly affect SDE2 splicing and expression level in heterozygous carriers compared with control animals. Initiation of translation at the closest in-frame methionine codon would truncate the SDE2 precursor by 83 amino acids, including the cleavage site necessary for its activation. Finally, no homozygote for the G allele was observed in a large population of nearly 29,000 individuals genotyped for the mutation. The low frequency (1.3%) of the derived allele in the French population and the availability of a diagnostic test on the Illumina EuroG10K SNP chip routinely used for genomic evaluation will enable rapid and efficient selection against this deleterious mutation.

Identification of copy number variation in French dairy and beef breeds using next-generation sequencing.

BACKGROUND: Copy number variations (CNV) are known to play a major role in genetic variability and disease pathogenesis in several species including cattle. In this study, we report the identification and characterization of CNV in eight French beef and dairy breeds using whole-genome sequence data from 200 animals. Bioinformatics analyses to search for CNV were carried out using four different but complementary tools and we validated a subset of the CNV by both in silico and experimental approaches. RESULTS: We report the identification and localization of 4178 putative deletion-only, duplication-only and CNV regions, which cover 6% of the bovine autosomal genome; they were validated by two in silico approaches and/or experimentally validated using array-based comparative genomic hybridization and single nucleotide polymorphism genotyping arrays. The size of these variants ranged from 334 bp to 7.7 Mb, with an average size of ~ 54 kb. Of these 4178 variants, 3940 were deletions, 67 were duplications and 171 corresponded to both deletions and duplications, which were defined as potential CNV regions. Gene content analysis revealed that, among these variants, 1100 deletions and duplications encompassed 1803 known genes, which affect a wide spectrum of molecular functions, and 1095 overlapped with known QTL regions. CONCLUSIONS: Our study is a large-scale survey of CNV in eight French dairy and beef breeds. These CNV will be useful to study the link between genetic variability and economically important traits, and to improve our knowledge on the genomic architecture of cattle.

Genetic variability of the activity of bidirectional promoters: a pilot study in bovine muscle.

Bidirectional promoters are regulatory regions co-regulating the expression of two neighbouring genes organized in a head-to-head orientation. In recent years, these regulatory regions have been studied in many organisms; however, no investigation to date has been done to analyse the genetic variation of the activity of this type of promoter regions. In our study, we conducted an investigation to first identify bidirectional promoters sharing genes expressed in bovine Longissimus thoracis and then to find genetic variants affecting the activity of some of these bidirectional promoters. Combining bovine gene information and expression data obtained using RNA-Seq, we identified 120 putative bidirectional promoters active in bovine muscle. We experimentally validated in vitro 16 of these bidirectional promoters. Finally, using gene expression and whole-genome genotyping data, we explored the variability of the activity in muscle of the identified bidirectional promoters and discovered genetic variants affecting their activity. We found that the expression level of 77 genes is correlated with the activity of 12 bidirectional promoters. We also identified 57 single nucleotide polymorphisms associated with the activity of 5 bidirectional promoters. To our knowledge, our study is the first analysis in any species of the genetic variability of the activity of bidirectional promoters.

Flt3 ligand improves the innate response to respiratory syncytial virus and limits lung disease upon RSV reexposure in neonate mice.

Respiratory syncytial virus (RSV) causes severe bronchiolitis in infants worldwide. The immunological factors responsible for RSV susceptibility in infants are poorly understood. Here, we used the BALB/c mouse model of neonatal RSV infection to study the mechanisms leading to severe disease upon reexposure to the virus when adults. Two major deficiencies in neonatal lung innate responses were found: a poor DCs mobilization, and a weak engagement of the IFNI pathway. The administration of Flt3 ligand (Flt3-L), a growth factor that stimulates the proliferation of hematopoietic cells, to neonates before RSV-infection, resulted in increased lung DC number, and reconditioned the IFNI pathway upon RSV neonatal infection. Besides, neonates treated with Flt3-L were protected against exacerbated airway disease upon adult reexposure to RSV. This was associated with a reorientation of RSV-specific responses toward Th1-mediated immunity. Thus, the poor lung DCs and IFNI responses to RSV in neonates may be partly responsible for the deleterious long-term consequences revealed upon adult reexposure to RSV, which could be prevented by Flt3-L treatment. These results open new perspectives for developing neonatal immuno-modulating strategies to reduce the burden of bronchiolitis.

Identification of large intergenic non-coding RNAs in bovine muscle using next-generation transcriptomic sequencing.

BACKGROUND: The advent of large-scale gene expression technologies has helped to reveal in eukaryotic cells, the existence of thousands of non-coding transcripts, whose function and significance remain mostly poorly understood. Among these non-coding transcripts, long non-coding RNAs (lncRNAs) are the least well-studied but are emerging as key regulators of diverse cellular processes. In the present study, we performed a survey in bovine Longissimus thoraci of lincRNAs (long intergenic non-coding RNAs not overlapping protein-coding transcripts). To our knowledge, this represents the first such study in bovine muscle. RESULTS: To identify lincRNAs, we used paired-end RNA sequencing (RNA-Seq) to explore the transcriptomes of Longissimus thoraci from nine Limousin bull calves. Approximately 14-45 million paired-end reads were obtained per library. A total of 30,548 different transcripts were identified. Using a computational pipeline, we defined a stringent set of 584 different lincRNAs with 418 lincRNAs found in all nine muscle samples. Bovine lincRNAs share characteristics seen in their mammalian counterparts: relatively short transcript and gene lengths, low exon number and significantly lower expression, compared to protein-encoding genes. As for the first time, our study identified lincRNAs from nine different samples from the same tissue, it is possible to analyse the inter-individual variability of the gene expression level of the identified lincRNAs. Interestingly, there was a significant difference when we compared the expression variation of the 418 lincRNAs with the 10,775 known selected protein-encoding genes found in all muscle samples. In addition, we found 2,083 pairs of lincRNA/protein-encoding genes showing a highly significant correlated expression. Fourteen lincRNAs were selected and 13 were validated by RT-PCR. Some of the lincRNAs expressed in muscle are located within quantitative trait loci for meat quality traits. CONCLUSIONS: Our study provides a glimpse into the lincRNA content of bovine muscle and will facilitate future experimental studies to unravel the function of these molecules. It may prove useful to elucidate their effect on mechanisms underlying the genetic variability of meat quality traits. This catalog will complement the list of lincRNAs already discovered in cattle and therefore will help to better annotate the bovine genome.

Construction and validation of a novel dual reporter vector for studying mammalian bidirectional promoters.

Regulation of gene expression plays important role in cellular functions. With the development of sequencing techniques, more and more genomes are available and genome-wide analyses of genomic structures that may affect gene expression regulation are now possible. Analyses of several genomes have found a class of regulatory regions that contain elements that initiate transcription of two different genes positioned with a head-to-head arrangement in two opposite directions. These regulatory regions are known as bidirectional promoters. Although bidirectional promoters have been known for years, recent genome-scale studies have shown that the regulation of the expression of up to 10% of the genes are controlled by bidirectional promoters. These findings are based mostly on computational work and only a limited number of putative bidirectional promoters have been experimentally validated. Developing methods to study bidirectional promoters will allow researchers to understand how these regions are regulated and the roles that divergent transcription plays in the expression of genes. Here, we have developed a novel dual-fluorescence reporter gene vector to study the transcriptional output of mammalian bidirectional promoters. We demonstrate that this vector is capable of expressing reporter genes under the control of bidirectional promoters, using the known human OSGEP/APEX bidirectional promoter.

Bovine TWINKLE and mitochondrial ribosomal protein L43 genes are regulated by an evolutionary conserved bidirectional promoter.

TWINKLE is a mitochondrial DNA helicase playing an important role in mitochondrial DNA replication. In human, mutations in this gene cause progressive external ophtalmoplegia and mitochondrial DNA depletion syndrome-7. TWINKLE is well conserved among multicellular eukaryotes and is believed to be a key regulator of mitochondrial DNA copy number in mammals. Despite its involvement in several diseases and its important function in mitochondrial DNA metabolism, nothing is known about the regulation of the expression of TWINKLE. We have analysed the 5'-flanking genomic region of the bovine TWINKLE gene and found it was localised adjacent to the MRPL43 gene in a head-to-head orientation, suggesting that both genes are regulated by a shared bidirectional promoter. The bovine 75-bp long intergenic region shows substantial homology across different species and contains several conserved putative transcription factor binding sites. A TATA box, however, was lacking. Using a dual fluorescent reporter system and transient transfection assays, we have analysed the bovine intergenic region between TWINKLE and MRPL43. This small genomic fragment showed a bidirectional promoter activity. As the TWINKLE/MRPL43 bidirectional promoter tested was highly conserved, it is likely that the results we obtained here in cattle may be extended to the other species.

Transcriptomic analysis of host immune and cell death responses associated with the influenza A virus PB1-F2 protein.

Airway inflammation plays a major role in the pathogenesis of influenza viruses and can lead to a fatal outcome. One of the challenging objectives in the field of influenza research is the identification of the molecular bases associated to the immunopathological disorders developed during infection. While its precise function in the virus cycle is still unclear, the viral protein PB1-F2 is proposed to exert a deleterious activity within the infected host. Using an engineered recombinant virus unable to express PB1-F2 and its wild-type homolog, we analyzed and compared the pathogenicity and host response developed by the two viruses in a mouse model. We confirmed that the deletion of PB1-F2 renders the virus less virulent. The global transcriptomic analyses of the infected lungs revealed a potent impact of PB1-F2 on the response developed by the host. Thus, after two days post-infection, PB1-F2 invalidation severely decreased the number of genes activated by the host. PB1-F2 expression induced an increase in the number and level of expression of activated genes linked to cell death, inflammatory response and neutrophil chemotaxis. When generating interactive gene networks specific to PB1-F2, we identified IFN-γ as a central regulator of PB1-F2-regulated genes. The enhanced cell death of airway-recruited leukocytes was evidenced using an apoptosis assay, confirming the pro-apoptotic properties of PB1-F2. Using a NF-kB luciferase adenoviral vector, we were able to quantify in vivo the implication of NF-kB in the inflammation mediated by the influenza virus infection; we found that PB1-F2 expression intensifies the NF-kB activity. Finally, we quantified the neutrophil recruitment within the airways, and showed that this type of leukocyte is more abundant during the infection of the wild-type virus. Collectively, these data demonstrate that PB1-F2 strongly influences the early host response during IAV infection and provides new insights into the mechanisms by which PB1-F2 mediates virulence.

Insights into gene expression profiling of natural resistance to coccidiosis in contrasting chicken lines.

Coccidiosis is a parasitic disease with major economic impact, one of whose main causative agents is Eimeria tenella. Chicken breeds display variable natural resistance to this disease. Unravelling the genetic bases of such variations could provide new clues for protection strategies. Transcriptomic experiments were conducted comparing resistant (Fayoumi) and susceptible (Leghorn) lines. Caecum and caecal tonsils were analysed. A global increase in differential gene expression following infection was observed for caecum comparisons, whereas a global decrease following infection was observed for caecal tonsils.Gene lists for infected tissues display 40 genes in common across breeds, 20 of which were specific to infected tissues. Among these specific genes, 9 belong to the 100 more differentially expressed genes of the infected caecum comparison. Gene expression networks were constructed in parallel, identifying highly connected genes. Comparing information from differential gene lists and gene network analysis allows one to highlight potential pivotal genes in the infection process, one of which was located in a putative significant QTL region for infection associated lesions.

Evolution of major milk proteins in Mus musculus and Mus spretus mouse species: a genoproteomic analysis.

BACKGROUND: Due to their high level of genotypic and phenotypic variability, Mus spretus strains were introduced in laboratories to investigate the genetic determinism of complex phenotypes including quantitative trait loci. Mus spretus diverged from Mus musculus around 2.5 million years ago and exhibits on average a single nucleotide polymorphism (SNP) in every 100 base pairs when compared with any of the classical laboratory strains. A genoproteomic approach was used to assess polymorphism of the major milk proteins between SEG/Pas and C57BL/6J, two inbred strains of mice representative of Mus spretus and Mus musculus species, respectively. RESULTS: The milk protein concentration was dramatically reduced in the SEG/Pas strain by comparison with the C57BL/6J strain (34 ± 9 g/L vs. 125 ± 12 g/L, respectively). Nine major proteins were identified in both milks using RP-HPLC, bi-dimensional electrophoresis and MALDI-Tof mass spectrometry. Two caseins (ß and αs1) and the whey acidic protein (WAP), showed distinct chromatographic and electrophoresis behaviours. These differences were partly explained by the occurrence of amino acid substitutions and splicing variants revealed by cDNA sequencing. A total of 34 SNPs were identified in the coding and 3'untranslated regions of the SEG/Pas Csn1s1 (11), Csn2 (7) and Wap (8) genes. In addition, a 3 nucleotide deletion leading to the loss of a serine residue at position 93 was found in the SEG/Pas Wap gene. CONCLUSION: SNP frequencies found in three milk protein-encoding genes between Mus spretus and Mus musculus is twice the values previously reported at the whole genome level. However, the protein structure and post-translational modifications seem not to be affected by SNPs characterized in our study. Splicing mechanisms (cryptic splice site usage, exon skipping, error-prone junction sequence), already identified in casein genes from other species, likely explain the existence of multiple αs1-casein isoforms both in SEG/Pas and C57BL/6J strains. Finally, we propose a possible mechanism by which the hallmark tandem duplication of a 18-nt exon (14 copies) may have occurred in the mouse genome.